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PROTON TISSUE DOSE FOR THE BLOOD FORMING ORGAN 
IN HUMAN GEOMETRY: ISOTROPIC RADIATION 

By G. S. Khandelwal* and John W. Wilson 
Langley Research Center 

SUMMARY 

A computer program is described which calculates doses averaged within five 
major segments of the blood forming organ (BFO) in the human body taking into accoimt 
seK- shielding of the detailed body geometry and nuclear star effects for proton radiation 
of arbitrary energy spectrum (energy less than 1 GeV) and isotropic angular distribution. 
The dose calculation includes the first term of an asymptotic series expansion of trans- 
port theory which is known to converge rapidly for most points in the human body. The 
result is always a conservative estimate of dose and is given as physical dose (rad) and 
dose equivalent (rem). 


INTRODUCTION 

A necessary requirement for radiation shield and dosimeter design and for space 
mission analysis is a reliable method of calculating the anticipated dose distribution in 
the human body for the pertinent radiation environment. Customarily, the human body is 
approximated by a simple geometric object (such as a sphere, slab, or cylinder) with 
resultant disagreement among the various approximations (ref, 1) and consequent dis- 
agreement on shield and dosimeter design and impact on mission objectives. 

Described herein is a computer program which calculates proton dose averaged 
over five major segments (upper limbs, lower limbs, upper trunk, lower trunk, and skull) 
of the blood forming organ (BFO) by treating the human body geometry in detail (refs. 2 
and 3) but assuming isotropicity of the incident primary particles. The calculation uses 
an approximate form of transport theory in which nuclear star effects are incorporated 
(refs, 4 and 5). The output of the program is in terms of physical dose (rad) and dose 
equivalent (rem). 


^Professor of Physics, Old Dominion University, Norfolk, Va. 


SYMBOLS 


d(x) 


fx(t) 

R(t,E) 

t 


X 

<f>{n,E) 


dose at the point X, rad or rem 


incident proton energy, MeV 


areal density distribution function, 


g/cm2 


dose response at depth t for protons of energy E, 

proton/cm<5 

depth from surface or areal density, g/cm2 

vector locating point within the body, g/cm2 

vector denoting direction of motion of radiation, dimensionless 

proton differential fluence spectrum, protons/cm^-MeV-sr 

THEORY 


An asymptotic expansion for the dose at a point X is derived in reference 5. The 
first term of the expansion is an accurate approximation of the dose given by 

R|t(n),E]0(n,E) dfidE (1) 

where the errors are second order in the expansion parameter 6, which is the beam 
width. In equation (1), 0 (o,e) is the proton fluence spectrum in the direction with 
energy E, R(t,E) is the fluence-to-dose conversion factor at the depth t for normal 
incident protons on a slab, and t^n) is the distance from the body surface to the point X 
along the direction Q. When the proton fluence spectrum 0 (o,e) is isotropic, equa- 
tion (1) can be written as 

d(x) = cj){E) R|t(s2),E] ds! dE (2) 

since 0 is no longer dependent on ^that is, 0^0, = 0(E)^. 

It is advantageous to replace the integral over by an integral over the depth t 
of tissue. This can be accomplished by the correspondence between the fractional solid 
angle and the depth t of tissue. Thus, the above integral can be written as 
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( 3 ) 


d(x) = 47T r 0(E) r R(t,E) fx(t) dt dE 

•^0 


with 


fx(t) dt = 1 (4) 

where fx(t) is the areal density distribution function for the point X. The quan- 
tity fx(t) dt is the fraction of the solid angle for which the distance to the surface 
from the point X lies between t and t + dt. The areal density distribution func- 
tions of Billings et al, (table 1) are discussed in reference 3 (see also refs, therein 
and in particular their ref. 23) and differ markedly from their original data in refer- 
ence 2. The function R(t,E) is given in references 4 and 5, In the program, the 
integrals over areal density and energy are performed using the Gauss-Legendre quad- 
rature formula (ref. 6). 


PROGRAM DETAILS 

The data for the program are self-contained except that the user supplies a sub- 
routine to generate the primary fluence spectrum (protons/cm2-MeV-sr), A sample 
spectrum is supplied with the program for initial checkout. The output for this sample 
calculation is given in tables 1 and 2. 

In the program, an intermediate function Rx(E) is calculated by evaluating 

Rx(E) = r R(t,E) fx(t) dt (5) 

^0 

after which the dose is evaluated as 
/^\ pEmax 

d(^x) - 47T Rx(E) 0(E) dE (6) 

where Ej^^x restricted to 1 GeV. A brief explanation of the principal subroutines is 
given below. 

There are three basic functions to be generated for the calculation. These are the 
areal density distribution fvmction, the fluence-to-dose conversion factors for normal 
incidence on a slab, and the incident fluence spectrum. Two numerical integrations are 
then performed to evaluate first equation (5) and then equation (6). These calculations 
are performed in the following subroutines: 
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Subroutine 


Purpose 


APDF Generates areal density distribution for five segments of BFO and 

BFO average 

ARES Generates the fluence-to-dose conversion factors R(t,E) 

SPECT0R Generates the incident fluence spectrum (protons/cm^-MeV-sr) 

R0FE Generates the function Rx(E) given by equation (5). The entry 

point RTAB generates a table of Rx(E) for interpolation. 

BD0S Generates the dose d(x) at point X by computing the final inte- 

gral over energy as indicated by equation (6). 

The program requires no input other than the incidence spectrum supplied by the user, 

A listing of the program is found in the appendix. The language is FORTRAN IV 
and written for the CDC 6000 series computers. Although it was originally programed 
on a 60-bit machine, round-off error on shorter word machines should not be a problem 
for these calculations. The incident spectrum in the list is 



with a sample output from the program given in tables 1 and 2. This sample calculation 
requires 112 seconds of CPU time on a CDC 6600 and 25000 (octal) locations. 

Lai^ley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 23, 1974. 
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APPENDIX 


PROGRAM LISTING OF THE DOSE CALCULATION 


PROG^^AM GSK ( INPUT ^OUTPUT ^ TAPE3 = OUTPUT ) 

COMi>/ION/ALPHA/Al.P 
real C05E(6»2) 

ALP= 1 ® 

CALL initial 
CALL BL)OS(UOSE) 

PRINT 1 

PRINT 2 ^ (OOSE ( I ) ^ I = 1 « 6 ) 

PRINT 3 W DOSE ( 1 ) ^ 1=7 ^ 12) 

1 FORMAT (30X<^ ARMS LEGS HI TRUNK LO TRUNK 

1 LL AVE BFO*// > 

2 FORMAT ( 30X^ --- REM-^^- 96 E 10« 1 // ) 

3 FORmAT(3ux^* RAD* ^ 6 i:i 1 0* 1 //// ) 

FNP 


SUBROUTINE INITIAL 
REAL AD (2H ) .PDF (6 ) 

data ad /I » 1 .1 ft 32 5 . 1 « 5 US . 1 ® 02 . 2 ® 3 . P ® . 3 « 3 . A « . 4 * 8 . E « 7E . 6 ® 9E 

1 1 Oft OS* 12ft. I4»E.l 7*^5»2l » *2S»“^*30»'r::.3S»E«i44ft ^S3® .F)3«5.7S« .91 » 
A = 1 ft 

call intapof 

PRINT 4 

4 FORMAT ( IHl .///// > 

PRINT S 

5 FORMAT!* THE BILLINGS AND LANGLEY ARRAl DENSITY D I STR I HL 

1 NOTIONS ARE* *//// ) 

PRINT 3 

3 FORMAT!* THICKNESS ARMS LEGS 

INK LO rPUrNX SKljLL AvE iiFO*.//) 

DO 1 I =1 .25 
CALL APDF ( Af) ! I ) .PDF ) 

1 PRINT 2 . AD ! I ). .PDF 

2 FORMAT !7E 15ft 3) 

PRINT 4 
CALL RTAB!A) 
return 
END 


s<u 


m8«35 . 
5/ 

'TION FU 
HI TRU 



u u 


APPENDIX - Continued 


SUBROUTINE BOOS (DOSE) 

C THIS SUBROUTINE COMPUTES THE DOSE AT POINT X BY COMPUTING THE 

C integral over energy fTHE INTEGRAND BEING THE PRODUCT OF 

C THE INCIDENT FLUENCE SPECTRUM AND THE DOSE RESPONCE 

C AT POINT X. RSUBX OF E 

REAL ANS ( 6f 2 ) ^ A (6 ^ 2 ) ^DOSF ( 6 ^ 2 ) 
external POX 
THK= 1 ® 1 01 

call am I ( thk #FM I N ) 

EMIN=EMIN+«5 
E0-TH = EM I N+10 , 

Ml =M2=23 

CALL iMGAUSS(EMIN,E0TH.M1 ,D0SE^P0X^ A, 12) 

CALL MGAUSS (FOTH,1000®,M2,ANS^POX^A,12) 

DO 1 L= 1 f 6 
DO 1 I =1 f 2 

1 DOSE (L^ I ) =D0SE (L , I )+ANS (Lt I ) 

RETURN 

END 


SUBROUT I NE POX ( X , F > 

THIS SUBROUTINE COMPUTES THE INTEGRAND NEEDED IN 
COMPUTING OOSEE D(X) AT POINT X 
REAL F ( 12 ) ^R ( 1 2 ) 

E = X 

call spec TOR (SPECTRA U:. ) 

SPfctRA = SPECTRA^'<' 1 2 • F66 
call ROFL(R.E) 

DO 1 I = 1 » 1 2 
1 F ( I ) =R ( I ) *SPECTRA 
RETURN 
END 


SURROUT I NE SPECTOR ( SPECTRA ^ E ) 

C THp SUBROUTINE GENERATES THE INCIDENT FLUENCE 

C SPECTRUM (P/CM2-“SR-MFV ) 

C O M M O N / A L P H A / A I . P 
bPECTRA=l «/E^^ALP 
RETURN 
END 
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APPENDIX - Continued 


SUBROU TINE ROFE ( R i T ) 

c THIS Subroutine generates rad ano pem oose responses at point 

C R SUB X OF F FOR FIVE SFGf’-^iENTS EFO AND f^FO AVERAGE. 

C THF entry point RTAR GENERATES A TAHL^^ RESPONSES AT 

C POINT X» R SUE X OF Efl FOR I NTE RPOL A T I ON 

COMMON /ENERGY/ XENER 

REAL R ( 12 ) ^ RT 1 2 ) ( 4 .SO -» IE ) ^ TAB ( EC ) 

REAL ANS ( 1 2 ) » A ( 1 2 ) 

DA TA T Ah/34 » ,40 * , 5'j ® , 60 ® , 7 , riu « , 90 ® ^ 1 ' ' C » , 1 1 ’ • ^ 1 20 « , 1 30 « , 1 BO * < 1 70 
1,2 O'O « , 2 bi : e , 3 o’ N * 9 4 ■ ’ r; » ,5 0 , 6 0 « , 7 0 ® , B ' : ’ • , 9^'"^ O © , 1 00 • / 

external fox 

DO 1 (’1 1 = 1,12 

1 01 call spl I n ( tab ,RT ( 1 , I ) ,NPT , C ( 1 , 1 , I ) , T ,R ( I ) ) 

RETURN 
.ENTRY RTAB 
NPT=2J 

DO 1 I = 1 9 NPT 
XENER=TAE ( I ) 

ENFR= XENER 
UENER=1 010© 

CALL RANGE ( R « ENER , UFNER , A , B ) 

TMIN=1 « 101 
FOTH=R 
C0TH=R 
TOTH=R 

IF (FOTH©GE©43©837 )FOTH = 4 3© 837 
IF(T0TH©GF©5® )T0TH=5© 

M 1 =?0 

call MGAUSS ( TM I N , TOTH, Ml , AMS , FOX , A , 1 2 ) 

no 11 J= 1 » 1 2 

11 RT ( I , J ) =ANS ( J ) 

IF (TOTH-b© ) 1 ,3, 1 

3 IF (COTH©GE©20© )COTH = 20. 

CALL MGAUSS ( TOTH , COTH , M 1 , ANS , FOX , A , 12 ) 

DO 12 J=1 , 12 

12 RT ( I , J ) =RT ( I , J l+ANS ( j ) 

IF (COTH-20© ) 1 , 4 , 1 

4 CALL MGAUSS(COTH,FOTH,vii ,ANS,F0X, A, 12 ) 

DO 13 J = 1 , 1 2 

13 RT ( I , J ) =RT ( I , J ) +AMS ( J > 

1 CONTINUE 

DO 1 T's J = 1 , 12 

1 0 CALL SPLCN (TAd,RT(l ,J),NPT,C(1^1 ,J) ) 

RETURN 

END 
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APPENDIX - Continued 


SUf^ROUTINE FOX(X,F) 

THIS SUBROUTINE COMPUTES THE INTEGRAND FOP PAO AND REM DOSE 
RESPONSES at point X, R SUB X OF E 
REAL F (6), PDF (6) 

COMV.ON /ENERGY/ XEMER 
T = X 

ENFR=XENER 

CALL ares ( T ^ ENER 9 RFiY ,RAD ) 
call APDF(T,P0F) 
no I I = 1 9 6 
F ( 1+6 ) =RAD*PDF ( I ) 

1 F ( I ) =REM*PDF ( I ) 

RETURN 

END 


SU3R0UT I ME APOF ( AD 9 PDF ) 

C THIS PROGRAM GENERATES THE BILLINGS AND LANGLFY AREAL DENSITY 
C distribution FUNCTION PDF ( I ) FOR FIVE SEGMENTS OF l3FO (1=195) AND 
C THE AVFRAGE OVFR BFO (1=6) 

REAL ADT (E5 ) iPQFT (6^35)9 TPDF (259 6 ) 9 0(49^596) 9 PDF ( 6 ) 9 FRAC ( 6 ) 
pATA ADT/1 « 1 9 1 .325,1 ,5/5,1 . G 2 9 2 « 3 , 2 « 75 , 3 . 3 , 4 . 9 4 , O 9 5 « 75 9 6 , 95 9 B . 35 9 
1 1 J ® .) 5 9 1 2 » 9 1 4 « 9 i ( .5921 ® 9 2 . 5 9 3 ® 5 9 3 6 « 5 9 4 4 . , 5 3 » 9 6 3 » 5 9 7 6 « 991 . 5 / 

DATA TP L^F/ «v)9/69 «/2B4 9 »56L ’ ^9 1.37^309 »B2^2 9 .9^9^, .57459 .5247 9 


1 

a 

7562 9 

2« 

1 11 6 9 

1 

« 

3v 75 

9 

9 80 0 3 9 

» 

4 1349 


, 1 5 72 

9 

9 

1 679 9 


« 1 5 0 0 9 


9 2 u 2 0 9 

2 

« 

1 4B0 , 


1 1 V 7 9 


9 

356 7 

, 

« ’ 4 4 5 9 

« 

0 1 74 9 


*0 124 

9 

a 

■JU7 1 9 

0 

« ' J 0 0 0 9 

0 

* coon. 

30 

a 

0 0 ^> 3 9 

V ' » 

' J \J sy , 


« 

B9 32 

9 

* B < 2 9 


9 I B / 9 

1 

9 6 0 8 9 

9 

2 9 

1 3B9 9 

2 

*07039 

2 

*6675, 

42 

9 

62 169 

2 9 

022A , 

1 

9 

72 4 7' 

9 1 

«470Y, 

« 

93 70 9 


e B 3 0 3 

9 

» 

4695 9 


® 2 5 6 1 9 


« 1 667 9 

6 

« 

14 94 , 

« 

06C'/9 


» 

0 3 74 

9 

9 ;'24B , 


U . 3 1 B 9 

C 

« 0 0 • 0 

9 

« 

0937, 


*2437 9 


*52779 

6 

« 

654 4 9 

e 

90 1 3 9 


9 

5971 

9 1 

* 1 B5 1 , 

I « 

15439 

1 

*3451 

9 

1 9 

23199 

1 

*4042 9 

1 

*28169 

71 

0 

ESS5 9 

9 

S* /IB, 


9 

959 3 

9 

9 6 294 9 

« 

53 1 7 9 


.3994 

9 

« 

2244 9 


*1144, 


« 0638, 

B 

e 

0350^ 9 

« 

029 1 9 


• 

1 2B 

9 

« 0 0 0 u' 9 


’ J w' C '.J 9 

O' 

» ' - 0‘ 0 

9 

» 

J 0 0 0 9 

0 

9 

0 

c5 

0 

0 

U 

, 0000, 

90 

9 

^’0 ^0 9 

9 

3 1 3 9 


9 

0 70 3 

9 

®2B1 39 

« 

70 45, 


9 92 3 1. 

9 

9 

E>4 589 


*8355, 


,96259 

1 1 

« 

'' ■■8^ )6 9 

« 

70 2 1 , 


9 

5984 

9 

9 2775 9 

, 

1 6 C 0 9 


9 1 295 

9 

, 

0703 9 


*0375, 


9 0 1 3 1 9 

1 

« 

C’336 9 

O' 9 

COCO , 


9 

0 000 

9 ' 

e ' 1 ^ , 

« 

0 1 37 9 


*3431 

9 

9 

74 71 , 


*7570 9 


*5/779 

0 

9 

3 5 8 6 9 

e 

4 6 B 2 9 


, 

4 3 1 0 

9 

95302 9 

9 

/ 7' 2 5 , 


9520 1 

9 

9 

33079 


*32 1 4 9 


*22589 

3 

e 

CEB 5 9 

e 

C" 1 75 9 


9 

4 3 

9 

9')1 37 9 

« 

1 3 9 9 


® OOBO 

9 

9 

0 0 43, 


*00579 

0 

« oooo. 

4 

a 

1 9 1 3 9 

» 

9 72 2 9 

I 

, 

1 'j22 

9 E< 

9 2 6 3 3 9 

-3 9 

8 1 9 8 9 

3 9 2 2 2 5 

9 

3 9 

97 25, 

4 

9 2 4 68 9 

4 

• 92 1 0 9 

57 

9 

14 6^9 

6 * 

786 7 9 


9 

7 72 6 

9 5 

»0B14- 

3® 

8 92 0 , 

3 

«4752 

9 

2 » 

5375 9 

1 

*8904 9 

1 

« 0986 9 

5 

« 

6050 9 

9 

4636 . 


9 

2532 


,1.3629 

s 

^■836 9 


9 0 60 9 

9 

0 ® 

0 000/ 






data FRAC/ «137^!-9 « 32969 6 O ' 69 el920t ®08l69l»0014/ 

DO 1 1=196 

call SPLIN( A0T,TPDF( 1 9 I ) , 25 9 C { 1 9 1 9 I) , AD 9 PDF ( I ) ) 

I F ( PO" { I ) . L T a 0 , ) PDF ( I ) = c . 

1 continue 
RF^TURN 

ENTRY INTAPDF 
DO 3 1=196 

DO 3 J= 1 925 

3 TPDF ( J 9 I ) nTPDF ( J 9 I )/ ( 1 00, ^FRAC ( I ) ) 

DO 2 1=196 

2 CALL SPLCN( ADT,TPDF( 1 9 I ) ,259 C( 1 9 1 9 I ) ) 

RETURN 

END 
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APPENDIX - Continued 


SUBROUT INE RANGE. ( R ^ ENER ajENER^ A « B ) 

C THIS SUBROUTINE COf^PUTES PROTON RAnGEi IN TISSUE 

9 IE (FNER-1 Oe ) 1 « 1 9 3 
3 IF (ENER-30 9 >6.69 1 1 

11 IF (FNER-1 CO* ) 1 3f i 39 1 2 

12 IF(ENER-3C0« >22922,23 
23 IF (FNER-bOO» >17917,25 
2b IF (ENER-UENER >69,69, 19 

1 A=:«002245 

B=1 ,698 
GO TO 60 
6 A=*0CM956 
B= 1 0 8 
GO TO 60 

13 A=»00 19775 
B= 1 .8 

GO TO 60 
22 A=*QU276624 
R=1 *729 
GO TO 60 
17 A=. 0054847 

R=1 *61 
GO TO 60 
69 A=*0 125462 
B=1 ,477 
GO TO 60 
60 CONTINUE 

r=A^FNER**B 

return 

19 WRITE (3,26 > FNER 

26 FORl^Af ( IHC, 1E15®6 > 

END 
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APPENDIX — Continued 


SUBROUTINF AMI(Fx,ENFR) 

C THIS SUBROUTINE CQivPUTES THE CUTOFF ENERGY FOR A GIVEN THICKNESS 

C OF TISSUE 

9 IF (EX-0 121 ) 1 0 1 ,3 
3 IF (EX-.89163 )6o 6o 1 1 

11 IF (EX-7.8726 ) 1 3o 13, 12 

12 I F (EX-b3,C68 )22 ,22,23 
23 IF (EX-1 21 .48 ) 1 7, 1 7, 25 

25 IF (EX-340, )69, 69, 5 

1 A=. 002245 

B= 1 ® 69B 
GO TO BO 
6 A=«001956 
B= 1 

GO TO 60 

13 

B=1 *8 
GO TO 60 . 

A= *00276604 
B=1 *729 
GO TO 60 
17 A=*00b4647 
B=1 *61 
GO TO 60 
69 A=*0l2b462 
B=1 *477‘ 

GO TO 60 

60 ( FX/A ) ( 1 »/R ) 

RETURN 
5 CALL EXIT 
END 

SUBROUT INE AEDUC ( A « F3 « L X U- NERP , LNLR ) 

THIS SUBROUTINE CO'^PUTES THF RESIDUAL ENERGY OF PROTON 
after IT TRAVERSES TISSUE OF A GIVEN THICKNESS 
AR = B 
aa=ex/a 

aar = aa^^--«- ( I »/AR ) 

2B IF (ENER-AAR)26^26927 

26 ENERP-0* 

RETURN 

27 ENFRP= ( ENFR^'^f ( AR ) -“A A ) ( 1 9 /AR ) 

RETURN 
END 
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APPENDIX - Continued 


SUf-'"ROUT I ME ANTFM (K.MtR . P0C3 > PTV.'O « POMQ , PO < P TV.' , PON . « ) 

THIS C3Ur.-F<OUT I NH GK N" RATt.TS THG VALUES OF THF PARAiVFTLRS 
OF- TKiL ANALYTIC FITS OF THF YONTE CARLO RESULTS 
D I MENS ION E ( 20 ) 9 Ao (20 ) , A2 ( 20 ) , AAO ( 2 0 ) ^ A A2 ( 2 0 ) *» C ( 4 ^ 2C ) 

DATA t / 30 « 4 O © ^ HO 9 6 ' ^ 70 » « B J e % 1 O O 9 , 1 2 ^ , 1 F O » , 1 70 « ^ 2 O 0 ® 9 

1 2 5 U 8 ^ 3 0 O ® * 4 , b o • s 8 60 u » ^7'3 0' » ? B C.! -O » % 9 n {.* » , 1 0 0 0 « / 

DATA AO / I « 4 7 « 1 » 4 7H » 1 » 4 0>0 © 1 ® ASF « I ® 49 © 1 » 4 9S , 1 * SO 9 1 © HOB t 1 ® 52 5 s 1 « 53 7 » 1 
1 ® 554 9I ®6il ®68<»1 e'^7^2«05^29 4 3,3©21 ^ 3 9 70 9 4 * 60 © 5 ® 75/ 

DA T A A-«2/u « O’ <9 « '-.jO ! ^ & ;ro2 ^ ® or; 3 9 © 0045 ^ « <'''CE> ^©01 <j » 0 1 9 © « 04 ^ © 058 © © 0729 1 « 057 
1 ^ »(>4- -S 9 *0305 © * 026 9 « 0225 -> • C' 1 77 , * 0 1 6 9 * o l 4 , * 0 I 2/ 

DATA A A 0/1 9 0 4? , 1 « 055 , 1 • 0.69 , \ •OBO * 1 .OQO* 1 • 1 , 1 . 1 2 * 1 . 1 48 , 1. 1 59 , 1 . 1 72 * 
1 1 ® 1 9 • < 9 1 » 2 C."^ 5 » 1 9 22 © I » 2 A © 1 « 2 7 8 ? i ® . -5 22©i®4?i ®45H?1® 56 ^ 1 ©67/ 

L'* A T A A .A 2 / « ] 9 © ® 1 / it'i t ® 1 6 « » 1 4 u © © I 4 3 © ® I 3 2 © 9 09 96© « 0 7 © ©0568© « O 4 9 5 © © O 4 3 B » « 

1 -'■375© *0314© » u 2 2 8 © » 0 1 El 8 © © 0 1 5 6 © ® ] 3 8 5 © « C.» 1 3 © * C 1 2 5 © « 0 1 18/ 

call SPLCN(E © AC ©29©C ) 

CALL SPL I N ( F © AO © 20 © C ©FNEP © POO ) 

CALL SPLCN( E ©A2 ©20 ©C) 

call SPLlNdi© A2©2 >’©C ©LNER ©PTuU) 

CALL SPLCN (E © AA J © 2^ © C ) 

call SPL I N ( F © AAo © 20 © C © ENER © PO ) 

call SPLCN(F©AA2©20©C ) 

call SPL I N ( F © AA2 © 20 © C © ENFR , PTU ) 

func=ptuo*r 

PGNQ= (EXP (FUNC ) ~POO ) /R 
FUNC“PT*w3fR 

PON= (EXP (FUNC )~P0 )/R 

RETURN 

END 


SI.J5R0UT I NFO AF ( STOPP © CALF © AO © bO ) 

THIS SUfTROUTINF COMPuT^^S THF QUALITY FACTOR AS A FUNCTION OF 

linear FiNE:RGY transfer 

IF (STOPP-35* ) 1 1 © 1 1 © 1 2 

QAL*^= 1 6 

rRETURN 

QALF = A0-^ (STOPP^* (F3Q ) ) 

RETURN 

FNO 


SURROUT INE A TOPP ( E NE R © S TO<^P ) 


THIS SUbROUTiNE: COMPUTES THE STOPPING PO';/ER FOR PROTON: IN TISSUE 

ZETA=ENER/938«21 1 

BLTAS= ( (ZETA* ( ZlTA-!-2 ® ) ) / ( ( ZETA4-1 « ) **2 ) ) 

V;6E= < ( 1 ® C2220 1 *1 0 (6 ) ) -JFlET AS ) / ( 1 * -BE T A S ) 

FBET = ALOG ( 'VF3E )-l3ETAS 

STOPP= « 30726 1 48->(- ( -2 9 237B342+ » 529726*Fi:;ET ) /BETAS 

RETURN 

END 



APPENDIX - Continued 


•SUBK'OUT I NE AROB (PROr", ^ F NEPP ^ UENKR ^ ENEP ? OFNEP ) 

C THIS SUL5ROUTINE CQfVPUT ES THE PRCoABILiTY OF A NUCLEAR REACTION NOT 

C OCCUR I NG BEFORE A GIVEN THICKNESS OF T I SSliF IS TRAVERSED 

ener=enerp 
I PRO= 1 

9 IF (ENER-l «6 )63«S4%64 

54 IF (ENER“-9«U )56^66<»6^- 

55 IF (ENER-BE® )67,67«58 

68 IF ( F“NFR “* 1 4 O « ) 6 9 ^ 6 « 70 

70 I F ( EMER-490 « ) 7 1 ^ 7 1 ^ 7F 
7?. I F (^NER-UFNFP ) 73® 73 ^ 19 
63 AO=t® 

AON=C ® 

AT'*'=0 ® 

ATH=Ee 
afo=e « 

GO TO 4 

66 A0= *95500944 
AON=«54524166 
ATW~^2 ® 3204599 
4TH=4* 1216447 
AF0=~2® 5963624 
GO TO 4 

67 A0=®B8724FF4 
AON= *9798068 
ATS'=1 3 •50256 1 
ATH=-1 75*26228 

AF0 = 4 18 *861 44 
GO TO 4 

69 A0=*2980451 
A0N=133*932 

Arw=-* 12453671 -«-l 0**^ (5 ) 

ATH= *55578849^1 U* ( 6 ) 

AFO = -* 96458432-^1 (7 ) 

GO TO 4 

71 AO=- *67860829 

AON=* 7705571 4* 1 0* (3 ) 

ATW = -* 1 7b05698-Jf lOe** (6 ) 

ATH= • 1 96 1 1 718-^ 1 0 ® ( 8 ) 

AF0 = ’-« 864 1 3965*1 0 •** (9 ) 

GO TO 4 

73 AO= *4793539 
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APPENDIX - Continued 


AON=.~» 1455bB36*10®** (4 ) 

ATW= » 1 4942072*1 ( 7 ) 

ATH=~e b56376 1 9*1 0«** (9 ) 

AFO= « 7578764 1 * 1 0 * ** ( 1 1 ) 

GO TO 4 

4 CFE= 1 ® ZENER 
CSE=CFE*CFE 
CTE=CSE*CFE 
CFO=CSE*CSE 

yf=ao+aon*cfe+atw*cse+ath*cte+afo*cfo 

IF ( IPRO-l )597»5 
7 IPR0=2 
YNU= YF 

ener=dener 

go to 9 

19 WRITE (3^26 ) ENFR 

26 FOR.VAT ( IHO^ 1 El 5®6 ) 

CALL EXIT 

5 PR0B=YF/YNU 
RETURN 

END 


SUBROUT I NE ARES (X INT ,LNER ^ RESF ^ v.'RESF ) 

THIS PROGRA^^ GENERATE'S THE FLUENCE TO DOSE CONVERSION FACTOR 
R(T.E) AS A FUNCTION OF TISSUE THICKNESS AND ENERGY 
EX=XiNT 
AQ=*06 
BQ=»B 

denfr^ener 
uener=ioio* 

56 continue 

call RANGE(R»ENLR^UENER« A,B ) 

call ANTLR ( ENER ^ F-OQ » F^T WO » P ONO ^ PO ^PT'a/^PGN 9 R ) 
call AEDUC (A ,EoEX«ENFRP»-NER ) 

IF (ENERP ) 34 . 33 9 34 

33 CONTINUE 
RES<^=0e 
WRFSF=0* 

enfr=cener 

RETURN 

34 CONTINUE 

CALL AROb (PROB ♦ENERP fUENER ♦ ENER ♦ DENER ) 

ener=enerp 

CALL AT0PP(ENER^ST0PP ) 

IF (STOPP)33»33» 1 

1 I F ( STOPP-2700 * ) 2 ♦ 2 ♦ 33 

2 CALL AF ( STOPP ♦ QALF ♦ AO ♦BO ) 

22 pes=prob*stopp*oalf 

COREO= ( P004-P0N0*X INT )/ (EXP (PTwO*X I NT ) ) 

C0REC= ( PO+PON*X I NT ) / (EXP (PTW*X INT)) 

RESF = PES*COREQ* 1 ®6E-“8 

IVRfsF = STCPP*PROB*COREC* 1 «6 E-B 

ENFR=DENER 

RETURN 

END 



APPENDIX - Continued 


SURROUT INE SPL I N ( X ^ Y , M , C ^ X I NT « Y I NT ) 

THIS SUBROUTINE INTERPULATES AMONG THE DATA POINTS TO PROVIDE VALVES 
ON SUB INTERVALS AFTER CONSTANTS ARE CALCULATED BY SUBROUTINE SPLCN 
DIMENSION X (BO ) 9 Y (BO ) ^ C (4f 50 J 
DATA IP/0/ 

IF (XINT-X (1 ) )7^ 1 «2 

1 YINT=Y ( 1 ) 

RETURN 

2 K=1 

3 IF (XINT-XCK+l ) )69 49B 

4 YINT=Y(K+1 ) 

RETURN 

5 K=K+1 

IF (M-K )7^7^ 3 

6 Y I NT = ( X (K4-I ) -X I NT )^ ( C ( 1 ^ :<)-«■( X (< + 1 ) - XI NT ) -?f-*2-f-C ( 3?K ) ) 
YINT=YINT+(XINT-X(K) )*(C(2^K)*(XINT-X(K) )**2 + C(4.K) ) 

RETURN 

7 PRINT 101 

101 FORMAT!* OUT OF RANGE FOR INTERPOLATION *) 

IP=|P+1 

IF(IP«GT*BO) GO TO 10 
RETURN 
10 CONTINUE 

A=ALOG ( “1 » ) 

END 
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APPENDIX - Continued 


SUBROUTINE 5PLCN(X.Y,M,C) 

THIS SUBROUTINE CALCULATES THE CONSTANTS TO BE USED IN A SPLINE 
FIT OF THE data 

D I MENS 1 ON X ( SO ) < Y ( bO ) • D ( bO ) » P ( bO ) » E ( bCJ > » C ( 4 • bO ) » A ( bO » J ) » B ( bO ) * Z ( bO 
1 ) 

MM=M— 1 

DO a K= 1 » MM 

D (l<) =X ( K+1 )-X (K ) 

P(K)=B (K )/6. 

2 E ( i< ) = ( Y (K + 1 ) -Y (K ) )/D (K ) 

DO 3 K=2,MM 

3 B (K ) =E (K ) -E ( K-1 ) 

A ( 1 ,2 ) =-l .-D ( 1 )/D (Z ) 

A ( 1 , 3 ) =D ( 1 )/D (2 ) 

A (2> 3 ) =P (2 )-P ( I ) +A ( 1 , 3 ') 

A ( 2»2 ) =2»* (P ( 1 )+P (^ ) ) -P ( 1 ) < 1 *2 ) 

A(2»3)=A(?,3)/A(2,2) 

3 (2 ) =B ( 2 )/A ( 2.2 ) 

DO 4 l< = 3.YM 

A(K,2)=2.+(P(K-1 )+P(K) )-P(K-1 )*A(K-1 .3) 

B(K)=B(K)-P(K-1 )+B(K-l ) 

A(K,3)=P(K)/A(K.2) 

4 B (K) =B (K ) /A (< ,2 ) 

Q = D(iV“2 )/D(M-1 ) 

A (M , 1 ) = 1 .+Q+A (M-2 . 3 ) 

A (M,2 )=-Q-A (M,l )*A(M-1 .3) 

B (M ) =B ( Y-2 )-A (M , 1 )*B (K-1 ) 

Z ( M ) =B ( M ) /A ( M ,2 ) 

MN=K-2 
DO 6 I = 1 . MN 
< = M- I 

6 2 (K ) =B <K )-A (K .3 )*Z(K+1 ) 

2 ( 1 )=-A ( 1 .2 (2 ) -A ( 1 , 3)*Z (3 ) 

DO 7 K=1 ,MM 
Q=1 »/(6.*D(K ) ) 

C ( 1 *K ) =Z (K) #Q 
C ( 2 .K ) =Z (K+1 ) *Q 

C ( 3. K ) = Y (K 5/D (K )-Z ('.< )*P (K ) 

7 C(4.K)=Y(K+1 )/D(K)-Z(K+l)*P(K) 

RETURN 

END 



APPENDIX - Concluded 

SUBWOUT INE: Jv/=GAUSS ( a « S'v.vi « FUNC ^ FOFX ^ NU-MBF-IR ) 

D I MENS ION U ( 5 ) R ( F ) ^ SUM ( 1 ) , FOFX { 1 ) 

DO 1 LL=1« NUMBER 

1 SUM{lL)=0»0 

IF ( A«FQ,B ) return 
U ( 1 )=*425B6PB30o09lB4 
U ( P ) = « 28330230F0353 76 
U ( 3 ) = • 1 6020F2 1 5H504H8 
U ( 4 ) = « 06746B3 1 65FSF08 
U ( F ) = « O 1 3 O 4 6 7 3 F 7 4 1 4 1 4 
R ( 1 ) =• 147762 1 1 23F7376 
R (2 ) =« 1 346333596^4 99B 
R ( 3 ) = « 1 0954 3 1 B 1 25 799 1 
R (4 ) =« 074725674575290 
R (5) =*033335672 154344 
F I NE=N 

DELTA=PrNE/(B^A ) 

DO 3 K= 1 *N 
XI -K-1 

F I NE=A+X I /DFLT A 
DO P 11=1.5 
UU = U( I I )/DELTA + F INE 
CALL FUNC (UU.FCrX) 

DO 2 JOYBOY=l .number 

2 SUM ( JOYEiOY ) =R ( I i ) -J^FOFX ( J0Yf30 Y ) +SUM ( JOYBOY ) 

DO 3 JJ = 1 .5 

UU= ( 1 e 0-U ( JJ ) )/DELTA + F I NF 
CALL FUNC (UU.FOFX) 

DO 3 NN=1 .NUMBER 

3 SUM (NN ) =R ( JJ ) -^t-FOFX < NN )+SUM ( NN ) 
no 7 T JK= 1 . MUMPER 

7 SUM ( I Jl< ) =SUM ( I JK ) /delta 
RETURN 
END 
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00 


THICKiMESS 


I^IOGE^OO 
i»3k:5E^00 
i®39!>E + 00 
i.92OE^00 
2«300E-^00 
2W30E + 00 
3.300E^Q0 
4«000E4-00 
4«800E^00 
3.750E+00 
6*950E'«'00 
8* 33UE-«-00 
iaOO^E+Ol 
i»200E^01 
l«450E'i-01 
i«750£+01 

zaooE^oi 

2«330E*i-01 

3.050E^01 

3.650E^0l 

4«400£^0I 

5*300Ei-01 

6^330h-5'01 

7.6a0E^'01 

9.130E^0l 


TABLE 1.- AREAL DENSITY DISTRIBUTION FUNCTIONS 


ARMS 

LEGS 

HI TRUNK 

7« 103E--03 

0. 

3.593E-03 

3*30iE^G2 

0. 

9.344E-03 

4«U31£-02 

0. 

2.023E-02 

9«993E-02 

2.710E-02 

2.509E-02 

6*033£-G2 

2.555E-C2 

3.456E-02 

6«91 1E^02 

2. 78 76-02 

2.289E-02 

4« I81E-02 

4. 881E-02 

4.525E-02 

3«8l9E-02 

6.489E-02 

4. 4266-02 

5.3046-02 

b.281E-02 

5.158E-02 

i.337E-ai 

8.093E-G2 

4.724E-02 

9.316E-02 

7.954E-02 

5.384E-02 

3. 8256-02 

6. 136E-C2 

4.914E-02 

3.009E-02 

3.231E-02 

6.091E-02 

1. 144E-02 

4.462E-02 

3.726E-02 

1.222E-02 

2.843E-02 

3.678E-02 

1.092E-02 

2.519E-02 

2.413E-02 

1.689E-02 

1.424E-02 

2.0396-02 

1.077E-02 

7.77CE-C3 

1.531E-02 

8. 7126-03 

3.058E-03 

8. 604E-03 

4. 127E-03 

4. 5 3 3E- 03 

4.387E-03 

3.239E-03 

1.842E-03 

2«446E-03 

1 » 2666—03 

1.1356-03 

1.377E-03 

9.0256-04 

7. 524E-04 

1.1166-03 

3. 167E-04 

9.648E-04 

4.908E-04 

0. 

0. 

0. 


LO TRUNK 

SKULL 

AVE BFO 

0. 

0. 

1.9106-03 

0. 

0. 

9.7086-03 

0. 

1.6796-03 

1. 1016-02 

0« 

4. 2056-02 

3. 2596-02 

0. 

9. 1566-02 

3.3156-02 

0. 

9.2776-02 

3.2186-02 

1 .6306-03 

7. 0806-02 

3. 9676-02 

3.661E-03 

4.3956-02 

4.2416-02 

1. 4656-02 

5.7386-02 

4.9146-02 

3.669E-02 

5.2826-02 

7.1376-02 

4.808E-02 

6.4986-02 

6.7776-02 

4. 9266-02 

8.8546-02 

5.7656-02 

4.3526-02 

6. 3746-02 

5.0746-02 

5.0136-02 

4.0536-02 

3. 8876-02 

5. 6756-02 

3.9396-02 

3.4706-02 

3.6576-02 

2.7676-02 

2. 5346-02 

3.1176-02 

7. 2066-03 

1.888E-02 

1.4456-02 

2.1456-03 

1.0976-02 

8.3336-03 

1. 7526-03 

6.8406-03 

6.7456-03 

1.6796-03 

4.6306-03 

3.6616-03 

1.7036-03 

2. 5286-03 

1. 9536-03 

9.8046-04 

1.3606-03 

6.8236-04 

5. 2706-04 

8.3486-04 

1.8756-04 

6. 9856-04 

6 .08 1 6-04 

0. 

0. 

0. 



TABLE 2.- DOSES FOR SAMPLE SPECTRUM 



ARMS 

LEGS 

HI TRUNK 

LG TRUNK 

SKULL 

AVE BFO 

REM 

7.5E--14 

24 4E-14 

3»2E-14 

4«3E-15 

4« IE-L4 

3490E-14 

RAD 

4.9E-14 

• 

m 

1 

2*0E-14 

3.3E-15 

2«3E-i4 

2«0E-14 
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